Short-time dynamic exponents of an Ising model with competing interactions 
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In this work the two-dimensional Ising model with nearest- and next-nearest-neighbor interac- 
tions is revisited. We obtain the dynamic critical exponents z and 9 from short-time Monte Carlo 
simulations. The dynamic critical exponent z was obtained from the time behavior of the ratio 
Fi = (M 2 ) mo =o/(M)^ x ~ i d//z , whereas the non-universal exponent 9 was estimated from the 

■ time correlation of the order parameter (M(0)M(t)} ~ t e , where M(t) is the order parameter at 
instant t, d is the dimension of the system and ((■ ■ ■)) is the average of the quantity (■ ■ ■) over 
different samples. We have also obtained the static critical exponents /3 and v by investigating the 
time behavior of the magnetization. 
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^ ■ The two-dimensional Ising model with nearest- and next-nearest- neighbor interactions J\ and J2, respectively, 
presents non-universal critical behavior for — J2/J1 > 1/2 jl], ^ |j, S, ^. The existence of a non-universal critical 
line for — J2/J1 > 1/2 was suggested for the first time by van Leeuwen in 1975 0, who investigated the fixed 
point structure and critical surface of two-dimensional Ising models. Subsequently, in 1977, Krinsky and Mukamel || 
conjectured that the critical behavior of the model should belong to the universality class of the two-component vector 
model. Also in 1977, Nightingale [j^] worked on finite size calculations of the model and showed that the exponents 
vary continuosly with the coupling ratio between nearest- and next-nearest-neighbor interactions. 
£j ■ The model is defined by the Hamiltonian 

where Si = ±1 are Ising spin variables and the notation and ((i, j)) denotes that each sum runs over nearest- 
and next-nearest-neighbors, respectively. At T — the ordering is ferromagnetic for — J2/J1 < 1/2, whereas it is 
C*"^ ■ super-antiferromagnetic (SAF) for — J2/J1 > 1/2 |10(| . In the SAF phase, a row (column) of up spins alternates with 
a row (column) of down spins, as depicted in Fig. jl]/ From both numerical j|[ [jl] and analytical |l(], |lj] approaches, 

• it was shown that in the plane of the reduced temperature kT/J\ against the competition parameter — J2/J1, the 
phase diagram presents ferromagnetic, paramagnetic and super-antiferromagnetic phases. The critical line between 

, the ferromagnetic and paramagnetic phases belongs to the universality class of the two-dimensional Ising model, 
whereas the critical exponents along the critical line separating the SAF ordered phase from the the paramagnetic 
phase are found to be non-universal. Static critical exponents on the SAF-Paramagnetic critical line were obtained 
by Monte Carlo simulations [l3| , high temperature series expansion || , renormalization group calculations [|], [yj , 
finite-size scaling |)|, cluster variation || and coherent-anomaly methods || fl5|| . It was observed that such critical 
exponents vary as a function of the coupling ratio —J2/J1, thus confirming non-universality in this model. 

Today it is well known that, far from equilibrium, the short-time relaxation of the order parameter follows a 
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universal scale form M(t) ~ mot |16 , where M(t) is the order parameter at instant t (measured in Monte Carlo 
Steps per Spin - MCS), too = M(Q) is a (small) initial order parameter value and 9 is a dynamic critical exponent 
associated to the anomalous increasing of the magnetization after the quenching of the system. Short-time Monte 
Carlo simulations of the model performed in a recent paper by Ye et. al ]l7| showed that the dynamic critical exponent 
9 is also non-universal along the SAF-Paramagnetic transition line, once it depends on the ratio —JijJ\. 

To our knowledge, for this model there are no previous attempts to obtain the dynamic critical exponent z, which 
is defined as r ~ £ z , where r and £ are time and spatial correlation lengths, respectively. Therefore, in this work we 
obtain the critical dynamic exponent z of the model from the time evolution of the ratio |L8| 

mo—0 j.d/z fn\ 

where M(t) is the order parameter at instant t, d is the dimension of the system and ((■ ■ •)) is the average of the 
quantity (• ■ •) over different samples with initial order parameter value toq. This ratio has proven to be useful in 
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determining the exponent z, according to recent results obtained for the 2D Ising model, q = 3 and q = 4 states Potts 
models |l8j an d at the tricritical point of the Blume-Capel model fl9|| . In this technique curves of ln(F 2 ) against 
ln(i) lay on the same straight line for different lattice sizes, without any re-scaling in time, resulting in more precise 
estimates for z. 

In this work we also reobtain the dynamic exponent 8 along the non-universal critical line using the time correlation 
of the order parameter |^(| 

(M(0)M(t))~t e , (3) 

and compare our results with those obtained in ref. | fl7|| , where the exponent 9 was obtained from the scale form 
M{t) ~ mot 8 . The advantage in the use of Eq. (||) is that one does not need to fix precisely the initial order parameter 
value Too- The only requirement is that (too) = 0. Contrarily, the scale form M{t) ~ m i e demands sharply prepared 
initial states with a precise value of too, besides a delicate limit too — » 0. 

We have also calculated the static critical exponents (3 and v of the model. Such exponents are obatined indirectly 
because, from numerical simulations, we are able to obtain only the ratios of exponents (3/vz, 1/vz and A = (d — 
2(3/u)/z, as discussed in Ref. In order to obtain the former ratio {f3/vz), we have used the scale relation |2l| 

M{t)~t-Pl v \ (4) 

according to which the order parameter mo = M(0) = 1 decays in the short-time evolution of the system. On the 
other hand, the ratio 1/vz was obtained from the scale relation 

9 K hiM(i, K )U =0 ~i 1/,/z , (5) 

which is obtained by differentiating the quantity In M (i, k) in relation to n at the critical point, where k = (T — T c ) /T c 
and T c is the critical temperature. 
From the scaling relation |2lj 

M 2 (t)~t x , (6) 

we have also calculated the exponent 
The exponent v is given by 
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whereas the exponent (3 is given by 



vz I \vz 



Finally, from the exponents z, v and A, the critical static exponent j3 may also be obtained from 

0=Z(d-z\). (10) 

The layout of this paper is as follows: In section |l| we define the order parameter of the model and we give details 
about the numerical simulations performed. In section |ni| we present the results obtained for both dynamic and static 
critical exponents. Finally, in section |y| we briefly discuss the main results of this work. 



II. SHORT-TIME MONTE CARLO SIMULATIONS 



Simulations were performed in two-dimensional lattices with periodic boundary conditions. The spin states were up- 
dated using one-spin-flip heat-bath algorithm. For each value of the competition parameter —J2/J1 > 1/2 considered, 
we used the corresponding critical temperature given in ref. @ . 



3 



FIG. 1: Example of one possible super-antiferromagnetic (SAF) ordering. In this example, rows of up spins alternate with 
rows of down spins. The other possible configuration of this ordered state corresponds to columns of up spins alternating with 
columns of down spins. 

Once the simulations were performed along the line between the SAF and Paramagnetic phases, the order param- 
eter M considered in this work is that corresponding to the SAF ordering. The SAF phase may be seen as two 
interpenetrating sublattices, each sublattice been formed by alternating rows, as the example depicted in Fig. [I]: One 
sublattice (say, sublattice A) corresponds to rows + + + + +, whereas the other sublattice (say, B) corresponds to 

rows . By inverting the signs of all spins in each sublattice, one obtains another configuration in which 

sublattice A is formed by down spins and sublattice B corresponds to up spins. Thus, this ground state is twofold. We 
call attention of the reader to the existence of another twofolded SAF ordering, which consists of columns of up spins 
alternating with columns of down spins. Thus, the SAF ground state is fourfold, what implies that the transition line 
SAF-Para is of kind 4— >1. Following Zittartz p2[ , this phase transition is candidate to exhibit non-universal behavior. 

The order parameter M is defined as M — (M A — AI B )/2, where 

N N 

mA = EE ( n ) 

x=l y=l 
N N 

x— 1 y— 1 

In order to deal with square sublattices, we used lattices of 2N = 144 rows and N = 72 columns. In the absence 
of a magnetic field, the symmetry M A = —M B is expected. Thus, when preparing initial states for numerical 
simulations we have taken this expected symmetry into account, i.e., we prepared initial t = configurations such 
that M A (0) = —M B (0). For the exponent z, one sees from Eq. (||) that two independent runs are necessary in order 
to obtain the ratio F^. One of them starts with mo = M(0) = 1 and the other run starts with mo = M(0) = 0. For 
mo = 1, we just set S x , y = 1 in sublattice A and S x ,y = — 1 in sublattice B. On the other hand, for too = 0, we 
set S x ,y — 1 in sublattice A and S XtV = — 1 in sublattice B, inverting randomly the sign of half of the spin variables 
in each sublattice. The ratios of exponents P/vz and 1/vz are obtained from Monte Carlo simulations with initially 
ordered states, i.e., initial states with too = 1. For the exponents 9 and A, we must have for the initial order parameter 
(mo) = 0. In order to obtain such an initial state, we give a probabitlity p = 0.5 for every spin in sublattice A to point 
up. Thus, we calculate M A (0) and impose M B (0) = — M A (0). The spin configuration in sublattice B is chosen in 
order to satisfy such imposition. Finally, the time interval used for the runs performed in this work was [1,150], and 
the final results for the dynamic critical exponents were extracted from a time interval [£1,^2] (1 < t± < t<i < 150) 
which maximizes the goodness of fit q. For example, suppose that from two time intervals and [£^, £2] one 

obtains q = 0.99 and q = 0.4, respectively. The critical exponent obtained from the time interval [£1,^2] is assumed 
to be correct. 

The error bars of the exponents 9, z, (3/vz and A, showed in the following section, are the standard deviation from 
the mean value obtained from five independent runs with different initial seeds for the random numbers generator. 
For the quantity (3/v and for the exponents (3 and v, which are obtained indirectly, the error bars were evaluated from 
usual error propagation. 
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• Time correlation of the order parameter 
XRef. [17] 
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FIG. 2: Non-universal behavior of the dynamic critical exponent 9, as a function of the ratio between nearest- and next-nearest- 
neighbor interactions — J2/J1. In the graph, the high of each X corresponds to the error bars found in Ref. jnj. The points 
with error bars shown explicitly are the results obtained in this work from the time correlation of the order parameter. 
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FIG. 3: Dynamic critical exponent z, as a function of the ratio between nearest- and next-nearest-neighbor interactions — J2/J1. 



III. RESULTS 



A. Dynamic critical exponents 

In this subsection we present the results obtained for the dynamic critical exponents 8 and z. In Fig. || it is shown 
the dynamic exponent 9 in function of the ratio — J2/J1, obtained from Eq. (^|). Non- universality is observed, since 
9 clearly depends on the ratio — J2/J1. This result is in accordance with that obtained in Ref. Jjl7||. In Fig. we 
present the dynamic exponent z as a function of —J2/J1, obtained from the ratio given in Eq. (|2|)!From Fig [3]onc 
can not affirm if the exponent z is universal or not, once the variation of z with the coupling parameter — J2/J1 is 
not so pronounced as observed for the exponent 9. Furthermore, since the major variation (or fluctuation) of z occurs 
close to the disorder point (— J2/J1 = 0.5), such fluctuation may result from a crossover effect. 

From both Fig. || and Fig. ||we observe that, as the coupling parameter —J2/J1 increases, the exponents 9 and z 
tend towards their known Ising values 0.19 • ■ • and 2.15 ■ • ■, respectively. 
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FIG. 4: Ratio P/vz, as a function of the ratio between nearest- and next-nearest-neighbor interactions — J2/J1. 
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FIG. 5: Ratio f5/v, as a function of the ratio between nearest- and next-nearest-neighbor interactions — J2/J1. 



B. Static critical exponents 

In this subsection we show, as a function of the ratio —J2/J1, the ratio of the exponents 0/vz, (3/v and A = 
(d — 2[3/v)/ z, and the static critical exponents j3 and v. In Figs. [|, || and pi as a function of the coupling parameter, 
we show the ratios (3/vz, (3/v and A. We observe that the behavior of these quantities as the coupling parameter 
— J2/J1 varies is quite similar to that observed for the exponent z. So, we are not able to affirm anything about 
non-universality of these quantities. Particularly, from Fig. ^| we observe that the ratio /3/u does not reach its known 
exact Ising value (3/v — 1/8 for large enough values of — J2I J\- We confirmed this result by Finite Size Scaling 
calculations combined with conformal invariance [^3) . 

The static critical exponent v, obtained from Eq. (||), is shown in Fig. ^ as a function of the coupling parameter 
—Jz/Ji- The non-universal behavior observed for this exponent was already found before using Monte Carlo |24[ ] and 
transfer- matrix j^] methods. In Fig. [7] we observe that the exponent v tends towards the known exact two-dimensional 
Ising value v = 1 as — Jij J\ increases. 

The non-universal behavior of the static critical exponent /3, evaluated from both Eqs. (^|) and (|To|), is shown in 
Fig. |^ as a function of the ratio — J2/J1. We observe the good agreement for the numerical values obtained for the 
exponent (3 from Eqns. (||) and In the second case, we note bigger error bars. This is due to the fact that the 
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FIG. 6: Exponent A = (d— 2f3/v)/ z, as a function of the ratio between nearest- and next-nearest-neighbor interactions —J%jJ\. 
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FIG. 7: Non-universal behavior of the static critical exponent v, as a function of the ratio between nearest- and next-nearest- 
neighbor interactions — J2/J1. 



evaluation of the exponent (3 from Eq. ( px| ) involves more exponents previusly obtained, and so the error propagation 
implies the bigger error bars observed. 
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FIG. 8: Non-universal behavior of the static critical exponent /3, obtained from Eqns. (0) (full circles) and ( p"o| ) (triangles up), 
as a function of the ratio between nearest- and next- nearest-neighbor interactions — J2/J1. 
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IV. SUMMARY AND DISCUSSION 

In this work we study the dynamic critical exponents of the two-dimensional Ising model with competing nearest- 
and next-nearest-neighbors interactions. We obtained the dynamic exponent z along the transition line between the 
SAF and Paramagnetic phases by using a method recently proposed based on the mixing of initial conditions |L8). 
To our knowledge, this work is the first attempt to obtain the exponent z for this model. We have also calculated the 
dynamic critical exponent 9 using the time correlation scale form of the order parameter [pof and our results are in 
good agreement with those obtained in ref . [ jl7[ . Non- universality was also observed for the static critical exponents 
(5 and v. In particular, as the ratio — J2/J1 increases, we observed that the dynamic critical exponents tend towards 
their known Ising values. However, the ratio fijv does not tends towards its known exact Ising value (3/v = 1/8 as the 
coupling ratio — J2/J1 increases. We have confirmed this result by means of Finite Size Scaling caculations combined 
with conformal invariance [E3f . 
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